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We address several different Casimir experiments where theory and experiment disagree. First 
out is the classical Casimir force measurement between two metal half spaces; here both in the 
form of the torsion pendulum experiment by Lamoreaux and in the form of the Casimir pressure 
measurement between a gold sphere and a gold plate as performed by Decca et al.; theory predicts 
a large negative thermal correction, absent in the high precision experiments. The third experiment 
is the measurement of the Casimir force between a metal plate and a laser irradiated semiconductor 
membrane as performed by Chen et al.; the change in force with laser intensity is larger than 
predicted by theory. The fourth experiment is the measurement of the Casimir force between an 
atom and a wall in the form of the measurement by Obrecht et al. of the change in oscillation 
frequency of a 87 Kb Bose-Einstein condensate trapped to a fused silica wall; the change is smaller 
than predicted by theory. We show that saturation effects can explain the discrepancies between 
theory and experiment observed in all these cases. 
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I. INTRODUCTION 

The Casimir force in the purest geometry is extremely 
fascinating: There is an attractive force between two par- 
allel, perfectly reflecting metal plates embedded in vac- 
uum; this geometry was treated by Casimir in his clas- 
sical paper [l| of 1948. One may say that the fact that 
Maxwell's equations have solutions in vacuum, that elec- 
tromagnetic waves need no medium to propagate in, gives 
rise to these forces. One way to calculate the forces is to 
calculate the energy in the system and then take minus 
the derivative with respect to separation, d, between the 
plates. At zero temperature the energy may be obtained 
as a summation of the zero-point energy of all the elec- 
tromagnetic normal modes of the system. The normal 
modes change when d changes and so do the zero-point 
energies. At finite temperature the actual population of 
the modes affects the forces. 

If the idealized perfectly reflecting metallic plates are 
replaced by plates of real materials the boundary con- 
ditions are different. There is still a force, but it is 
now different. For real metals the vacuum modes pene- 
trate to some extent into the plates and furthermore new 
modes appear, modes that are bound to the surfaces of 
the plates; these are evanescent waves, so-called surface 
modes 0|. For large separations the vacuum modes, or 
propagating waves, dominate and the force, at least at 
zero temperature, is identical to the force in the ideal- 
ized system. For smaller separations the surface modes 
dominate and the force acquires different characteristics; 
in this range the force is called the van der Waals force. 

In a pioneering work Sparnaay 0] tried to verify the ex- 
istence of the Casimir force experimentally between two 
parallel plates. However, the uncertainties were of the 
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same order of magnitude as the force itself so the exper- 
iment was non-conclusive. The interest in the Casimir 
force virtually exploded a decade ago. This increase in 
interest was spurred by the torsion pendulum experiment 
by Lamoreaux [4j|, which produced results with good 
enough accuracy for the comparison between theory and 
experiment to be feasible. This stimulated both theorists 

0. 0.0, US, II3| and experimentalists [U El El El and 
the field has grown constantly since then. Another rea- 
son for this increase of interest in Casimir forces is the 
huge shift of general focus of the science community into 
nanoscience and nanotechnology where these forces be- 
come very important. 

A setback came immediately. Theory and experiment 
agree quite well for low temperatures, but at room tem- 
perature, where most experiments are performed there 
are serious deviations. Each new type of experiment 
has lead to new puzzling discrepancies between theory 
and experiment. In the torsion pendulum experiment by 
Lamoreaux Q and in the micromachined torsion oscil- 
lator experiment by Decca et al. the Casimir force 
and Casimir pressure, respectively, were found to agree 
much better with the theoretical zero temperature results 
than with the room temperature results although both 
experiments were performed at room temperature. The 
problems appeared when one improved the treatment of 
the dielectric properties of the metals by inclusion of dis- 
sipation Q; this can be done either by using tabulated, 
experimental optical data or by using the Drude dielectric 
function. The origin of the dissipation is carrier scatter- 
ing against impurities and lattice imperfections; at finite 
temperature also carrier-phonon scattering contributes. 
Removing dissipation from the Drude dielectric function, 

1. e., using the so-called Plasma model has the effect that 
the discrepancies go away. One was facing the dilemma 
that improving the theory led to worse results. There 
were also claims that the theory did not obey the Nernst 
heat theorem, the second law of thermodynamics. This 
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could, however, be shown not to be the case 
There still remained a problem for perfect crystals where 
the dissipation comes from phonon scattering, only. This 
was, however resolved by inclusion of spatial dispersion 
[ljj; it turned out that spatial dispersion had a very small 
effect on the Casimir force, but it removed the remaining 
problem with the Nernst heat theorem for perfect crys- 
tals. 

The sphere and plate geometry used in these two ex- 
periments has the benefit compared to the two plate ge- 
ometry used in Refs. 0, EH that one avoids the align- 
ment problem. However, one has to rely on the so-called 
proximity force approximation (PFA) |2y] in the interpre- 
tation of the experiments and that could be a source of 
error. Fortunately, the PFA, although based on a rather 
shaky foundation, seems to work quite well (2ll |. 

To get a better insight into what goes wrong, or where, 
it would be very useful to have experiments performed in 
the whole temperature range from very low temperature 
up to room temperature or even higher, keeping every- 
thing else fixed. This is very difficult to do in practice. 
Chen et al. [HJ came up with a very clever idea about 
how to collect complementing information. By using a 
gold sphere and a laser irradiated semiconductor wafer 
one could study how the force changed with carrier den- 
sity in the wafer, keeping everything else fixed. 

Here one was facing another setback. The change in 
force in the presence of light compared to in the absence 
of light was larger than expected. In the previous ex- 
periments one could get agreement between experiment 
and theory by neglecting the dissipation. In this exper- 
iment it did not matter if dissipation was included or 
not. One could get good agreement if one completely 
neglected the very few carriers, present in the absence 
of light but kept the contribution from the laser excited 
carriers in the presence of light. Thus, it seemed to be 
quite different reasons for the discrepancies in the two 
types of experiment. 

Another experiment [13] where the Casimir force can 
be measured is based on Bose-Einstein condensates 
(BEC) near a surface. For a harmonically trapped BEC 
placed at a distance d from the surface of a substrate 
the centre-of-mass oscillation frequency changes due to 
the surface-atom potential. Experiments were performed 
[HI on a 87 Kb BEC trapped to a fused silica wall. Theory 
predicted a too big a change in frequency. Good agree- 
ment was obtained if one neglected the contribution to 
the dielectric properties of the fused silica from the very 
small amount of thermally excited carriers in the material 

m. 

Thus, in the last two experimental setups it seems as if 
the dielectric properties of the carriers in a semiconduct- 
ing system with very small carrier concentration have to 
be modified. This has been done in two recent publica- 
tions [IE H3 but if this resolves the problems is not so 
clear. 

In this work we put forward what we think is the so- 
lution to all the problems discussed above, or at least 



the first step towards a solution. It is applicable to all 
the discussed experimental setups. The idea is that the 
discrepancies are due to saturation effects. We have pre- 
viously discussed these ideas in a very short and con- 
densed article [2|| . Here, we make a much more detailed 
presentation. In Sec. Owe present the formalism that is 
common to all three experimental geometries. Then all 
three geometries are treated one by one, starting by two 
metal plates in Sec. IIIIl followed by the semiconductor 
and metal plates in Sec. llVl and ending by the atom wall 
geometry in Sec. [V] Finally, we end with summary and 
discussions in Sec. IVI1 



II. BASIC FORMALISM 

The formalism we use is based on the electromagnetic 
normal modes of the system Q. At zero temperature 
the force, F(d), between two objects as a function of 
separation, d, is given by the derivative of the energy, E: 



F{d) 



dE(d) 
dd ' 



(1) 



where the energy is the sum of the zero-point energies of 
the modes, 



(2) 



In a simple system with a small number of well-defined 
modes this summation can be performed directly. In the 
general case this is not so. Then we can rely on an exten- 
sion of the so-called Argument Principle [l,l22j to find the 
results. Let us study a region in the complex frequency 
plane. We have two functions defined in this region; one, 
<p(z), is analytic in the whole region; one, f(z), has poles 
and zeros inside the region. The following relation holds 
for an integration path around the region: 
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where zo and the zeros and poles, respectively, 

of function f{z). In the Argument Principle the func- 
tion <p(z) is replaced by unity and the right hand side 
then equals the number of zeros minus the number of 
poles of the function f{z) inside the integration path. If 
we choose the function f(z) to be the function in the 
defining equation for the normal modes of the system, 
/ (uji) — 0, the function ip (z) as Hz/2 , and let the con- 
tour enclose all the zeros and poles of the function f(z) 
then Eq. ([3]) produces the energy in Eq. (J2j) - By using 
this theorem we end up with integrating along a closed 
contour in the complex frequency plane. In most cases it 
is fruitful to choose the contour shown in Fig. [TJ We have 
the freedom to multiply the function f(z) with an arbi- 
trary constant without changing the result on the right 
hand side of Eq. ^ . If we choose the constant carefully 
we can make the contribution from the curved part of the 
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Figure 1: Integration contour in the complex z-plane suited 
for zero temperature calculations. Crosses and circles are 
poles and zeros, respectively, of the function f(z). The ra- 
dius of the circle is let to go to infinity. 



contour vanish and we are only left with an integration 
along the imaginary frequency axis: 
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Figure 2: Integration contour in the complex z-plane suited 
for finite temperature calculations. Crosses and circles are 
poles and zeros, respectively, of the function f(z). The small 
semi circles are centered at the poles of the coth function in 
the integrand. The radius of the large semi circle is let to go 
to infinity. 



E(r) = — I diulnf(iLu). 



(4) The coth function has poles on the imaginary z-axis and 
they should not be inside the contour. The poles are at 



where the result was obtained from an integration by 
parts. At finite temperatures the force is [30] 



f( r ) = --T-, 
clr 

where Helmholtz' free energy is 
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and all residues are the same, equal to 2/h(3. We inte- 
grate along the imaginary axis and deform the path along 
small semicircles around each pole. The integration path 
is illustrated in Fig. [2l The integration along the axis 
results in zero since the integrand is odd with respect to 
u>. The only surviving contributions are the ones from 
the small semicircles. The result is 



We make the observation that we may use the generalized 
Argument Principle but now with In [2 sinh ((3hz/2)]/(3 
instead of Hz/2 for <p(z) in the integrand. There is how- 
ever one complication. This new function has poles of its 
own in the complex frequency plane. We have to choose 
our contour so that it includes all poles and zeros of the 
function f(z) but excludes the poles of <p(z). The poles 
of function ip(z) all are on the imaginary frequency axis. 
We use the same contour as in Fig. [H but now let the 
straight part of the contour lie just to the right of, and 
infinitesimally close to, the imaginary axis. We have 
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Since the summand is even in n we can write this as 



3f = j J2 ' ln f ( iu} n);oJ n = n = 0, 1, 2, . . . , (10) 



where the prime on the summation sign indicates that 
the n — term should be multiplied by a factor of one 
half. This factor of one half is because there is only one 
term with \n\ = in the original summation but two for 
all other integers. When the temperature goes to zero 
the spacing between the discrete frequencies goes to zero 
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Figure 3: Integration contour in the complex z-plane. Crosses 
and circles are poles and zeros, respectively, of the function 
f(z). The contour should be extended to include the whole 
positive part of the real axis. This contour is useful if one 
wants to find the contributions from different mode types 



and the summation may be replaced by an integration: 

oo 
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and we regain the contribution to the internal energy 
from the interactions, the change in zero-point energy of 
the modes. 

In all the three geometries we discuss in this work the 
interaction energy per unit area, V (d) , can at zero tem- 
perature be written on the form [2] 

oo 

^ (<*) = £ E/ (12) 
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where d is the distance between the objects, k the two- 
dimensional wave vector in the plane of the plate(s), ft 
the area of a plate, and / (k, to) = is the condition for an 
electromagnetic normal mode in the particular geometry. 
The integration is along the imaginary frequency axis. 
At finite temperature the integration is replaced by a 
discrete summation over Matsubara frequencies, 
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Alternatively one may integrate along the real frequency 
axis, 

oo 

V = E Im / % ("0 + V2] In [/ (k, co)} , (14) 
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Figure 4: Dispersion curves for the modes between two gold 
plates in absence of dissipation. The frequencies are in units 
of lu s , the surface plasmon frequency. The solid straight line 
is the light dispersion curve in vacuum; the dashed (dotted) 
curves are TE (TM) propagating modes; the thin solid curves 
are evanescent TM modes; the thick solid curve is the lower 
boundary for transverse bulk modes in the plates. From Ref. 
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where n (lo) = [exp {h(3uS) — 1] 1 is the distribution func- 
tion for massless bosons. This form can also be used 
at zero temperature; then the distribution function van- 
ishes. These last results one arrives at by using the con- 
tour in the complex frequency plane shown in Fig. [3l The 
force per unit area, or pressure, is obtained as minus the 
derivative with respect to distance, F (d) = —dV (d)/dd. 

In all three geometries there are two groups of normal 
mode, transverse magnetic (TM) and transverse elec- 
tric (TE), each with a different mode condition func- 
tion. The interaction potential is a sum of two terms, 
V (d) = V™ (d) + V TE (d). In the two plate geometries 
the mode condition functions are 
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where the Fresnel amplitude reflection coefficients at an 
interface between medium i and j are 
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for TM modes (p-polarized waves) and 
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for TE modes (s-polarized waves) , respectively. We have 
let the objects be of medium 1 and 2 and let the sur- 
rounding vacuum be denoted by medium 0. The gamma 
functions are 



lj (k, w) = ^/fc 2 -e,H (cu/c) 2 ; j = 0, 1, 2. 



(18) 
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Figure 5: The thermal contribution to the Casimir energy 
from the four mode types between two gold plates. The results 
are also lumped together in various ways. All results are 
divided by the zero temperature Casimir energy for perfectly 
reflecting plates. Similar results have been obtained earlier 
for silver plates (3^ | 

Now we have all material we need for now. The inter- 
action potential for the atom wall geometry can also be 
derived from these relations using a limit procedure. We 
make that derivation in Sec. [V] 



III. TWO PARALLEL METAL PLATES 

The dispersion curves for the electromagnetic normal 
modes for two gold plates is shown in Fig. [4] . This 
figure is valid in neglect of dissipation in the plate ma- 
terials. The modes are propagating (evanescent) above 
and to the left (below and to the right) of the light dis- 
persion curve. The light dispersion curve is the straight 
diagonal line in the figure; it has slope unity with the 
chosen scaling of the axes. Note that there are no TE 
evanescent modes. When the system is allowed to have 
dissipation there are modes everywhere. Each original 
mode is replaced by a continuum of modes [3l| . Evanes- 
cent TE modes appear and the continuum extends all the 
way down to the momentum axis. These modes are the 
cause of all the problems with the thermal Casimir force 
in this geometry. 

In Fig. [5] we show the thermal correction to the 
Casimir energy. The contribution from all four mode 
types are given separately. These curves have been ob- 
tained from Eq. (fl4l) and the Drude dielectric function 
suited for gold has been used. We note that the TE 
evanescent waves give a negative contribution. For larger 
separations the total TE contribution saturates when it 
has completely eliminated the zero temperature TE con- 
tribution. 

In Fig. [6] we show the energy correction factor for the 
two-gold-plate geometry. The energy correction factor 
is the Casimir energy per unit area divided by the zero 
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Figure 6: Energy correction factor for two gold plates. The 
filled squares with error bars are the Lamoreaux' experimental 
Q values from the torsion pendulum experiment. The dashed 
curves are the perfect metal results. The thick solid curves 
are the results for real gold plates using Eq. Ill2[| for zero tem- 
perature and Eq. l]13[l for room temperature. The dielectric 
properties for gold was obtained from tabulated experimental 
optical data. The curves with circles are the results from our 
model calculations with different saturation parameters. 



temperature Casimir energy per unit area for two per- 
fectly reflecting metal plates, — fc7r 2 /(720<i 3 ) . The ex- 
perimental results 3 obtained by Lamoreaux are shown 
as squares with error bars. The larger the separation 
the bigger the errors. There is a cluster of experimental 
points near lfxm separation. Here the deviation between 
theory and experiment is over 20%. Theory and experi- 
ment are in clear disagreement. The proposed prescrip- 
tion has been to neglect the dissipation in the intraband 
part of the dielectric function but keep it in the interband 
part [33[. 

Let us now explain our view of what goes wrong in the 
theory of the thermal Casimir effect in presence of dissi- 
pation. The traditional theory relies fully on the concept 
of electromagnetic normal modes. These are assumed 
to be independent massless bosons. The possibility to 
excite one of these modes is assumed to be completely 
independent of how many modes are already excited. An 
excitation of a mode involves excitations of the charged 
particles in the system, electrons in the geometries stud- 
ied here. These are the sources of the fields. Now, the 
electrons are fermions and there is at most one electron 
in each particle state. An electron that is excited at one 
instant of time cannot immediately be excited again - 
the state is empty. The more modes that are excited 
the more difficult it is to excite new modes — there are 
saturation effects. It is a matter of logistics. In the theo- 
retical treatment this is not taken care of. In most cases 
this fact will not cause any problems, but sometimes it 
could. We think that the thermal Casimir effect is one 
such case. When dissipation is included each mode is 
replaced by a continuum of an infinite number (for an 
infinite system) of new modes. The distribution function 
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Figure 7: The thermal contribution to the Casimir energy 
from the TE evanescent modes between two gold plates for 
different damping parameters, circles. The corresponding to- 
tal corrections are shown as solid curves. Only the TE evanes- 
cent modes have been modified by the saturation. 



diverges towards zero frequency and the saturation ef- 
fects should appear here. This is very difficult to treat in 
a strict way. We use an approximation which is very easy 
to implement. We shift the distribution function in Eq. 
(fl4ll downwards in frequency, so that it never reaches the 
point of divergence, by adding a damping parameter, D, 



n (w) = [exp (hflw +D)-1]~ 



(19) 



The discrete frequency summation in Eq. (fT3| is the re- 
sult of the poles of the distribution function that all fall 
on the imaginary axis, see Ref. [2]. Our new distribution 
function has its poles shifted away from the axis the dis- 
tance D/h/3 into the left half plane. The new form of the 
interaction potential is 



k u> n 
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(20) 

Each term in the summation is replaced by an integral. 
For small D values it is enough to replace only the zero 
frequency term. 

The circles in Fig. |(6|) are the results with damping pa- 
rameters 0, 0.01, 0.1, and 1.0, respectively, counting from 
below. We have used Eq. (fl4|) with the modified distri- 
bution function to get the thermal correction. A word of 
caution is in place. This approach changes the popula- 
tion of modes at all energies, not just at the low frequency 
limit where saturation effects are predicted to appear. To 
limit the damage of this we have just taken into account 
the change in the contribution from TE evanescent waves. 
These contributions are the dominating ones at low ener- 
gies. In Fig. [7] we show the thermal correction from the 
TE evanescent waves, including the damping, as circles 
and the total correction as solid curves. 



Figure 8: Energy correction factor for two gold plates. The 
filled squares with error bars are the Lamoreaux' experimental 
[J] values from the torsion pendulum experiment. The dashed 
curves are the perfect metal results. The thick solid curves 
are the results for real gold plates using Eq. lfT2|) for zero tem- 
perature and Eq. l]13[l for room temperature. The dielectric 
properties for gold was obtained from tabulated experimental 
optical data. The curves with circles are the results from our 
model calculations with different saturation strengths. Here 
the distribution function for all mode types is limited to the 
value M. The Drude dielectric function was used in the cal- 
culation of the thermal correction 



An alternative way to treat the saturation is to simply 
introduce a cutoff in the distribution function. In Fig. 
[5] we show the results for the energy correction factor 
from using Eq. lfl4|) with a cutoff. The parameter M is 
the maximum allowed value of the distribution function. 
We have here limited the separation range to focus on 
the interesting cluster of experimental points near 1/xm 
separation. This cutoff method forces us to stay on the 
real frequency axis where the results are more sensitive 
to variations in the dielectric function and the numerical 
calculation more difficult to perform. 

The experimental result [HI] for the normalized Casi- 
mir pressure at 295 K is shown as dots in Fig. [9l The 
bars are the endpoints of the experimental error bars. 
The upper (lower) thick solid curve is the theoretical re- 
sult for zero temperature (295 K) calculated with Eqs. 
([T2"ll and (f]~3|) . respectively. The dielectric function on 
the imaginary frequency axis was derived from experi- 
mental tabulated optical data for gold. We note that the 
zero temperature result agrees much better with the ex- 
perimental result. The large negative thermal correction 
comes entirely from the TE evanescent modes (34|. All 
curves are normalized with the zero temperature Casimir 
pressure between two perfectly reflecting metal plates, 
heir 2 / (240 z 4 ^ . We have neglected surface roughness ef- 
fects. The circles are the results for different damping 
parameters from using Eq. lfl4|) with the modified dis- 
tribution function of Eq. (fl9|) in the contribution from 
the TE evanescent waves. To each set of circles corre- 
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Figure 9: Casimir pressure between two gold plates. The ex- 
perimental result from Ref. is shown as dots and the 
endpoints of the error bars are indicated by horizontal bars; 
the upper (lower) thick solid curve is the traditional theoreti- 
cal zero (room) temperature result; the circles are the present 
results, obtained from Eq. lfT4|l with the distribution function 
for TE evanescent waves modified according to Eq. fUjl , with 
damping parameters 0.01, 0.1, and 1.0, respectively, counting 
from below; the corresponding results obtained by shifting the 
zero frequency pole of the cosh function into the left of the 
complex frequency plane are shown as thin solid curves. 



sponds a thin solid curve. This curve is the result of 
using Eq. fl3j) . where just the zero frequency pole of the 
cosh function has been moved into the left half of the 
complex frequency plane and the corresponding term in 
the summation has been modified according to Eq. (|20| . 
We note that for the two lowest set of curves with small 
damping the two results agree. For very high damping 
there are deviations. These deviations have two reasons: 
One is that in the thin solid curves all mode types are 
affected by the damping; the second is that for strong 
enough damping more terms in the summation, leading 
to the sets of circles, should be modified. 
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Figure 10: The change in Casimir force, at 300 K, between 
a gold sphere and a silicon membrane with and without laser 
irradiation. The open squares with error bars are the experi- 
mental [13] result. The dashed curve with open circles is the 
theoretical result without saturation effects. The solid curve 
with filled (open) circles is our present result with D equal to 
0.01 (0.1). 



similar. Besides, it is now the TM modes that cause the 
problems. In the non-irradiated semiconductor material 
there is always some thermally excited carriers present. 
One postulated that for the non-irradiated semiconduc- 
tor one should completely omit the contribution to the 
dielectric function from the thermally excited carriers. 
That brought theory and experiment into agreement, see 
Fig. 10 in Ref. [12]. In the present Fig. [10] the solid 
curve with filled (open) circles is our saturation based 
result with D equal to 0.01 (0.1). We find that both 
these curves agree with the experiment within the exper- 
imental uncertainty. Here we have used Eq. (fl3| and 
just modified the zero frequency contribution according 
to Eq. |20l) . We have derived the dielectric function for 
the gold plate from experimental tabulated optical data. 
For the silicon with and without laser irradiation we have 
used the dielectric function given in Ref. [22]. 



IV. TWO PARALLEL PLATES, ONE METALLIC 
AND ONE SEMICONDUCTING 

The third experiment we consider here is the measure- 
ment of the Casimir force between a gold plate and a 
laser irradiated semiconductor membrane as performed 
by Chen et al. [§2]. They measured the change in force 
with the laser irradiation compared to without any irra- 
diation. The idea is to find out how the force varies with 
carrier concentration in the semiconducting membrane. 
The results are shown in Fig. [lOl The open squares with 
error bars are the experimental result. The dashed curve 
with open circles is the theoretical result for 300 K. The 
deviations are clear. In this geometry it is not enough to 
neglect dissipation to get agreement with experiment; the 
theoretical results with and without dissipation are very 



V. ATOM WALL GEOMETRY 

The Casimir force between an atom and a wall can be 
obtained from the results of the two plate geometry. One 
takes the limit when the thickness of one of the plates 
goes to zero and at the same time lets the material of 
the thin plate be diluted. We start with a half space of 
medium 1 and a layer of medium 2 emersed in vacuum 
(medium 0) separated by the distance d. The function in 
the condition for electromagnetic normal modes is then 
given by 



/=1 



-27od r io r Q2 + e 272rf2 ri r 2 o 
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(21) 



2 r 02^20 

where we have suppressed the arguments k and oj. We 
have furthermore divided the original function with the 
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one corresponding to infinite separation. This corre- 
sponds to choosing our reference system to be the one 
at infinite separation. There will be one function like 
this for TM modes and one for TE modes. Using that 
the thickness of layer 2 is small and that the material in 
this layer is diluted with a small reflection coefficient as 
a result leads to 



/ = 1 + e-^" d r w r 2 



1 



-272 d,2 



1 _ e -2y 2 d 2 ( rQ2 )2 
l + e - 2 7o<V loro2 2 72 d 2 . 



(22) 



We now let the dielectric function in the layer be £2 = 
1 + Anna (lu) wl| 4ira (w)/Qd,2. We have assumed that 
there is only one atom in the layer, with polarizability 
a(uj). At this point we need to specify the mode type. 
We start with TE modes: 



1 + e 



-2 7o d 7i -70 7o - 72 
7i + 7o 7o + 72 



2d 2l2 . (23) 



We furthermore have 



72 = J k 2 - (lu/c) 2 - 4ira (lu) (lu/c) 2 / Qd 2 

« 70 (k, lu) - 2na (lu) (w/cf / £ld 2 ~f (k,u). 

Thus, the result is 

jTE _ x + 1 c -2 70 d 7i - 70 2™ (m) (m/cf 
^ 7i + 7o 7o 



(24) 



(25) 



2d 2 72, (26) 



For the TM modes we have 

f™ w 1 + c - 2 ~t» d l1 ^ g i7o£27o -72 , 
71 + £ i7o £270 - 

and in the limit 

2na(Lu) hk 2 - (uu/cf 



TM 



1 , 1 c -2 70 d 7i~gi7o 

£1 7! + £170 



/ 

^ "I -I- - I -■(! "Ml 

(27) 

In this geometry we want the interaction energy, not the 
interaction energy per unit area, so we multiply with the 
area of the plate and find 



00 

1, J 



(28) 



+ nin/™(k,«j)] , 
for zero temperature and 

h 



V(d) 



(3Q 



["ln/ TB (fc,^„) 



k lj„ 



(29) 



+ filn/™(fe,^„)] 



for finite temperature, where 



n\nf TE (k,w) 



ae 



-2 7o rf7i -7o 2tt (cj/c)^ 



71 + 7o 



7o 



(30) 



and 



mnf IM (k,Lu) w ae 



-2 7o ri7l 



~U0 



2tv 



2k 2 - (lj/cY 



71 + £170 



7o 



(31) 

respectively. We note that the zero frequency contri- 
bution vanishes for the TE modes, but not for the TM 
modes. The force between the atom and the wall is 



F(d) 



ae ^2 70d 7i_7£ {uJn/c) 2 



71 +7o 



ae 



-2 70 d 7i -£i7o 
7i + £i7o 



2k 2 + (coJcY 



(32) 

where the first term is the TE contribution and the sec- 
ond the TM contribution. The suppressed arguments of 
the functions in the integrand are (k,iw n ). 

In the experiment by Obrecht et al. f24| one stud- 
ied indirectly the force between a rubidium atom and 
a dielectric substrate. This was done by measuring the 
collective oscillation frequency of the mechanical dipole 
mode of a BEC near enough to a dielectric substrate 
for the force to measurably distort the trapping poten- 
tial. In that work one studied both the equilibrium case 
when the substrate was of the same temperature as the 
surrounding environment and the non-equilibrium case 
when the substrate was of a different temperature. We 
concentrate on the equilibrium case when both the sur- 
roundings and the wall have the same temperature, 310 
K. The fractional change in trap frequency is defined as 
7a; = (loq — lu x )/luq in terms of the unperturbed trap fre- 
quency, luq, and ui x , the trap frequency perturbed by the 
force. This may be obtained as 



lx = 



-^ll*e (d)\, 



(33) 



where the expression for <I> e (d) is identical to Eq. lf32|) 
except that the integrand or summand has an ex- 
tra factor l\ (270a) g (2jqR x ). In the experiment a 
'i.oi')//;;; is the amplitude of oscillation, R x = 2.69/j.m is 
the Thomas-Fermi radius in the x-direction and m — 
1.443 x 10" 25 fc<7 is the mass of the rubidium atom. 
The function I\ (z) is the Bessel function and g (z) = 
15 [(3 + z 2 ) sinhz — 3zcoshz]/z 5 . These results were 
obtained in Ref. [HI after an analytical averaging proce- 
dure had been performed. 

In Fig. Ql] the experimental result [HI is shown as 
open squares with error bars. The upper (lower) curve is 
the theoretical result, without saturation, including (ne- 
glecting) the conductivity from the few thermal carriers 
in the silica wall. We see that also here the neglect of the 
contribution, to the dielectric function of the silica wall, 
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7 8 9 10 11 

d (flm) 

Figure 1 1 : Fractional change in trap frequency for a Rb atom 
near a silica wall versus separation in thermal equilibrium. 
The open squares are the experimental result (24l |. The upper 
(lower) curve is the theoretical result including (neglecting) 
the conductivity from the few thermal carriers in the silica 
wall. The circles are our present results for the D values 
1CT 10 , 1CP 11 , and 1CF 12 , respectively, counted from below. 

from the very few thermally excited carriers brings the 
theoretical result into agreement with experiment. This 
neglect is the postulated remedy in Ref. [25||. In this ge- 
ometry, just as in gold-plate silicon-wafer geometry, the 
TM modes cause the problems and it is not enough to ne- 
glect dissipation to get good agreement between theory 
and experiment. To include saturation effects we have 
just modified the zero frequency contribution in analogy 
with Eq. (|20ll . We note that in this experiment it is 
enough to have a damping parameter as small as 10~ 10 



to bring the theoretical result into agreement with exper- 
iment. We have assumed that the thermally excited car- 
riers in the wall material have the conductivity 100 s -1 
(~ 10~ 10 ohm~ 1 crn~ 1 '), which is the upper limit of the 
range given in Ref. [25||. Using smaller values leads to 
even weaker demands on the damping parameter. 



VI. SUMMARY AND DISCUSSIONS 

In summary we have proposed that saturation effects 
are responsible for the discrepancy between theory and 
experiment in several quite different Casimir geometries. 
It is quite difficult or even impossible to derive the sat- 
uration effects from first principles. We have treated 
saturation within a very simple calculation model and 
demonstrated that the problems may go away in all cases. 
The basic assumption is that the electromagnetic normal 
modes can not be fully treated as independent massless 
bosons. There are limitations on the number of excited 
modes and this limitation shows up at low frequencies 
where the distribution function for massless bosons di- 
verges in the case of non-zero temperature. We encourage 
experimentalists to design experiments where the occu- 
pation of these modes is directly investigated. 
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